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ABSTRACT 

Context. dSph galaxies are the least luminous, least massive galaxies known. Recently, the number of observed galaxies in this class 
has greatly increased thanks to large surveys. Determining their properties, such as mass, luminosity and metallicity, provides key 
information in our understanding of galaxy formation and evolution. 

Aims. Our aim is to provide as clean and as complete a sample as possible of red giant branch stars that are members of the Hercules 
dSph galaxy. With this sample we explore the velocity dispersion and the metallicity of the system. 

Methods. Stromgren photometry and multi-fibre spectroscopy are combined to provide information about the evolutionary state of 
the stars (via the Stromgren ci index) and their radial velocities. Based on this information we have selected a clean sample of red 
giant branch stars, and show that foreground contamination by Milky Way dwarf stars can greatly distort the results. 
Results. Our final sample consists of 28 red giant branch stars in the Hercules dSph galaxy. Based on these stars we find a mean 
photometric metallicity of -2.35 ± 0.31 dex which is consistent with previous studies. We find evidence for an abundance spread. 
Using those stars for which we have determined radial velocities we find a systemic velocity of 45.20 ± 1.09 kms~' with a dispersion 
of 3.72 kms -1 , this is lower than values found in the literature. Furthermore we identify the horizontal branch and estimate the mean 
magnitude of the horizontal branch of the Hercules dSph galaxy to be Vo = 21.17 ± 0.05, which corresponds to a distance of 147+, 
kpc. 

Conclusions. When studying sparsely populated and/or heavily foreground contaminated dSph galaxies it is necessary to include 
knowledge about the evolutionary stage of the stars. This can be done in several ways. Here we have explored the power of the C\ 
index in Stromgren photometry. This index is able to clearly identify red giant branch stars redder than the horizontal branch, enabling 
a separation of red giant branch dSph stars and foreground dwarf stars. Additionally, this index is also capable of correctly identifying 
both red and blue horizontal branch stars. We have shown that a proper cleaning of the sample results in a smaller value for the 
velocity dispersion of the system. This has implications for galaxy properties derived from such velocity dispersions. 

Key words. Galaxies:dwarf - Galaxies: fundamental parameters - Galaxies: individual: Hercules - Galaxies: kinematics and dy- 
namics - Galaxies: photometry 



1. Introduction 

In the past few years the number of known Milky Way 
satellites has increased considerably. Our Galaxy has gained 
at least ten newly recogni zed companions, a n d additional 
ones a wait confirmation (e.g. [Zucker et al.ll2 006a b: Walshetal] 
l2007t iBelokurov etafl 12009. |2007l 120081 I2009I) . These re- 
cently discovered satellites resemble the previously known 
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dwarf spheroidal g alaxies in many of their characteristics (see 
iGrebel et all 12003. for a detailed discussion of the properties 
of the classical dSph galaxies), but most of the new satellites 
are several magnitudes fainter than any dwarf galaxy known be- 
fore. Hence these objects are now often referred to as "ultra- 
faint dwarf spheroidal galaxies". Since most of the new dis- 
coveries were made using deep CCD sky surveys such as the 
Sloan Digital Sky Survey, which cover primarily the northern 
hemisphere, it seems highly likely that additional objects will be 
added once the southern sky is scann ed in a similar manner (e.g. 
iKeller et al.ll2007l:IWalsh et alJl2009h . 

One of the new discoveries is the H ercules dwarf spheroidal 
(dSph) galaxy (Be lokurov et alj 120071) . Hercules lies at a dis- 
tance of 132 ± 12 kpc from us, has an absolute V-band mag- 
nitude of about -6.6 ± 0.3, a V-band surface brightness of 
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only 27.2 + 0.6 mag arcsec~ 2 , and ap pears highly elongated 
dColeman et al.ll2007t iMartin et alj|2008l) . Its stellar mass is es- 
timated to be in the range of several xlO 4 M Q . The mass of the 
Hercules dSph galaxy, as inferred from line-of-sight radial ve- 
locity meas urements, is of the order of 10 7 M within the cen- 
tral 300 pc (Strigari et al.l2008h . This value is in good agreement 
with the seemingly ubiq uitous, common mass scale of the other 
Galactic satellites ( e.g. iGilmore et ail 120071: IWalker et ai]|2007l 
Strigari et al. 2008, and references therein). 

As would be expected from its low luminosity and low sur- 
face brightness, Hercules is a metal-poor dSph galaxy. From 
measurements o f the Ca n IR triplet lines in red giant stars, 
Simon & Geha d2007l) inferred a mean m etallicity of [Fe/H~| 
= -2. 27 ± 0.07 on the metallicity scale of ICarretta & Grattonl 
d 1997b. Using spectrum synthesis of Fe i lines. iKirbv et al.l 
(2008b) derived a mean metallicity of -2.58 ± 0.51 dex. Both 
studies found a wide range of metallicities among the red gi- 
ant stars in Hercules, confirming the trend known from brighter 
dSph galaxies that often exhibit spreads of 1 d ex and more (e.g. 
IShetrone et ai1l200lHKoch eT^l2006ll2007allbh . 

A detailed, high-resolution abundance analysis of two red 
giants in the Hercules dSph galaxy revealed that the enrichment 
in heavy elements proceeded inhomogeneously and that core- 
collapse supernova e were the primary contributors to the enrich- 
ment of Hercules dKoch et al.ll2008bl) . Evidence for such chem- 
ical inhomogeneities on small s cales has also been found in the 
more luminous dSph galaxies (Koch et alj|2008al) and in more 
massive dlrr galaxies dKniazev et al 

In our current study we explore the potential of Stromgren 
photometry for the study of the stellar content of ultra-faint dSph 
galaxies, such as the one in Hercules. In contrast to the usual 
broadband photometry, Stromgren photometry offers several 
advantages. It provides us with gravity-sensitive multi-colour 
indices useful for distinguishing different evolutionary stages 
including giant-dwarf discrimination. This is a very valuable 
option when trying to eliminate foreground dwarfs from a 
giant candidate sample to be used for subsequent spectroscopy, 
particularly when dealing with sparse, extended, ultra-faint 
dSph galaxies that tend to suffer from substantial Galactic 
foreground contamination. Moreover, Stromgren indices offer 
the possibility to estimate metallicities for red giant stars. This 
method has been considerably r e fined s ince the early calibration 
attempts by Grebel & Richtler (1992). As compared to other 
photometric estimates, it has the added ad vantage of prov iding 
age-independent metallicity estimates (e.g. Faria et al. 2003)- In 
terms of telescope time, intermediate-band photometry is a lot 
cheaper than spectroscopic surveys of faint giants. 

We have obtained both Stromgren photometry as well as 
spectroscopic observations of the Ca 11 IR triplet lines for stars 
in the field of the Hercules dSph galaxy. These observations en- 
able a full analysis both of evolutionary stage as well as radial 
velocities for these stars. It turns out that knowledge about the 
evolutionary stage of the stars is crucial for the construction of 
a clean sample of red giant branch (RGB) stars in the Hercules 
dSph galaxy. 

The paper is organized as follows: In Sect. [2] we describe 
the observations and data reductions for both the photometric 
and the spectroscopic observations. In Sect. [3] we present the 
colour-magnitude diagram as well as the measured radial veloc- 
ities. In Sect. [4] we show how the gravity sensitive Stromgren c\ 
index can be used to disentangle the Hercules dSph galaxy mem- 
bers from the foreground contamination, and Sect. [5] deals with 
membership determination based on radial velocities. Section|6] 
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Fig. 1. Position on the sky of our observations. The coordinates 
are given in decimal deg rees, epoch J2000. Central coordinates 
for the galaxy are from Mar tin et"aT1 d2008l) . The solid ellipse 
represents the core radius and the dotted ellipse the King pro- 
file Umiting radius of th e Hercules dSph galaxy as determined 
bv lColeman et al.l (120071) . Solid rectangles outline the four CCD 
chips in the WFC used for the photometric observations, x 
marks the fibre positions for the FLAMES observations. Here 
we show the fibres positioned on the stars only. The arrow in the 
bottom left hand corner indicates the direction to the centre of 
the Milky Way. 

summarizes how we define a Hercules member star. Section [7] 
provides a comparison with previous velocity determinations. In 
Sect. [8] we derive metallicities for the stars identified as mem- 
bers of the Hercules dSph galaxy using the Stromgren m\ index 
as well as from the measurements of the Ca 11 IR triplet lines, 
and compare the results with previous studies of the metallicity 
of the Hercules dSph galaxy. In Sect.|9]our results are discussed 
and they are summarized in Sect. [TO] 

2. Observations and data reductions 

In this section we detail the observations as well as the data re- 
duction for both photometric and spectroscopic observations. 

2.1. Photometry 

Intermediate-band Stromgren u, v, b and y photometry was ob- 
tained during six nights with the Wide Field Camera (WFC) on 
the 2.5-m Isaac Newton Telescope (INT) on La Palma. Of the six 
nights only three provided useful data due to bad weather con- 
ditions. The WFC consists of 4 2kx4k CCDs. The CCDs have 
a pixel size of 13.5 microns corresponding to 0.33" per pixel. 
Figure Q] shows the location and dimensions of the galaxy on the 
sky and the positions of the CCD chips of the WFC. 

The observations are summarised in Table Q] and all obser- 
yations were center ed at RA=16 /l 31 m 05 I and Dec=+12°47'18" 
(IMartin et al.l l2008). Typical seeing during the three good nights 
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Table 1. Summary of the photometric observations obtained 
with the Isaac Newton Telescope. 



Filter 


13 April 2007 
[min] 


14 April 2007 
[min] 


15 April 2007 
[min] 


Total 
[min] 


y 


1 x30 


1 x30 


1x30 


90 


b 


1 x30 


1 x30 


1 x30 


90 


V 


1 x30 


1 x30 


1 x30 


90 


u 


1 x30 


2x30 


2x30 


150 



Column 1 lists the filter. Column 2-4 list the number of 30 minutes 
exposures obtained for each filter during each of the three useful 
nights. Column 5 lists the total exposure time for each filter. 



Table 2. St andard stars from [Schuster & Nissenl dl988l) . except 
HD107853 dQlsenlll993l) . 



ID 


Hip 


V 


(b-y) 


mi 


C\ 


HD 100363 


56327 


8.648 


0.191 


0.139 


0.760 


HD 107853 




9.100 


0.321 


0.157 


0.472 


HD 108754 


60956 


9.006 


0.435 


0.217 


0.254 


HD 120467 


67487 


8.147 


0.728 


0.757 


0.088 


HD 134439 


74235 


9.058 


0.484 


0.224 


0.165 


HD 138648 


76203 


8.137 


0.504 


0.358 


0.290 


HD 149996 


81461 


8.495 


0.396 


0.164 


0.305 


HD 158226 


85378 


8.494 


0.386 


0.146 


0.316 


DM -05 3063 


51127 


9.734 


0.568 


0.461 


0.182 


DM -08 4501 


87062 


10.591 


0.452 


0.032 


0.274 


DM -12 2669 


43099 


10.230 


0.229 


0.094 


0.490 


DM -13 3834 


69232 


10.685 


0.415 


0.098 


0.183 


DM -14 4454 


81294 


10.332 


0.565 


0.469 


0.192 



Column 1 and 2 list the star ID and Hipparcos number, respectively. 
Column 3 lists the standard values adopted for the magnitude (note 
y = V) and column 4 to 6 the standard values adopted for the (b-y), 
mi and c\ indices, respectively. 



was about 1.3 arcsec. 

Multiple standard and extinction stars were observed each 
night. During the observations of the Hercules dSph galaxy, we 
observed 12 Stromgren standard stars chose n from the l ist in 
ISchuster & Nis sen (1988) plus one star from lOlsenl d 1993b . see 
Table [2] Two stars from the list of standard stars were observed 
several times during each night in order to sample the extinction 
for a large range of airmass. These stars will henceforth be re- 
ferred to as extinction stars. 

The observation of standard and extinction stars are used to find 
the zeropoint, extinction coefficients and colour terms, see Sect. 



2.1 .1 . Reduction of the photometric observations 

The images for the Hercules dSph galaxy and the standard stars 
were reduced with the Wide Field Survey Pi peline provided 
by th e Cambridge Astronomical Survey Unit (llrwin & Lewis! 
2001). The processing operations applied to the images were 
de-biasing, trimming, flatfielding, astrometry and correction for 
non-linearity. 
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Fig. 2. The residuals for the standard stars as a function of our 
final (v - y) colour. The residuals are given in the sense "our 
observed value" - "the standard value". Mean differences and <x 
as indicated. 



2.1 .2. Standard star photometry and establishing the 
photometric calibration 

We obtained aperture photometry for the standard and extinction 
stars using the task PHOT within the IRAfQ APPHOT package. 
The size of the aperture was determined individually for each 
star by plotting the measured flux as a function of increasing 
aperture size. The aperture at which the flux no longer increased 
was chosen as the aperture for that star (typically 4 - 5 x FWHM 
of the stellar psf). This curve-of-growth is used in order to 
maximise the signal-to-noise ratio (S/N) while measuring as 
much flux as possible from the star dHowellll 1989b . 

The measurements of the standard and extinction stars were 
used to establish the transformati ons needed to put our obser- 
vations on the standard s ystem of lOlsenl dl993i) . Note that the 
ISchuster & Nisseri) dl988l) stars are secondary standards in this 
system. See lOlsenl d 19951) for a discussion about the Olsen sys- 
tem as compared to the system est ablished bv lBondl d 19801) and 
lAnthonv-Twarog & Twarogl d 19941) . 

The first night out of the three useful nights did not give us re- 
liable photometry for the standard stars (i.e. there was no well 
defined linear relation for the transformation to the standard sys- 
tem for that night), and this night was thus excluded from the 
calibration. 

First, we derived preliminary extinction coefficients, and 
zeropoints, zu for each filter i by solving the following equation 
with a least-square fit 

m s = m + ki ■ X + Zi (1) 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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where m is the magnitude of the star and X is the airmass. The 
subscripts s, and i designates the standard magnitude, the ob- 
served magnitude and the filter, respectively. With preliminary 
zeropoints and extinction coefficients we then solved for the 
colour term, a,, using a least-square fit 



Table 4. Aperture corrections. 



= ra + at ■ (v - y) s + z/ 



(2) 



where z, is the residual from the linear fit. Note that here we use 
the data for both nights in order to make use of the colours for all 
standard stars. With preliminary colour terms we can now solve 
the full equation to obtain better estimates of Zi and k. 



= m + a t ■ (v - y) s + k" ■ X + z" 



(3) 



With these better estimates, Eq. (f2]) is again solved and Eq. 
© and OJ are then iterated until convergence in a, z and k is 
achieved. Typically convergence is reached within 6 to 7 itera- 
tions. The final zeropoint is given by z" in Eq. (O plus the resid- 
ual zeropoint from Eq. (f2]) from the last iteration. The extinction 
coefficients, colour terms, zeropoints and correlation coefficients 
between the uncertainty in z and k, p^, for each filter and night 
are listed in Table [3] We note that the uncertainty in zeropoint 
and extinction coefficient are strongly anti-correlated. 

The z, k and a obtained were used to tran sform our observa- 
tions onto the standard system of lOlsenl d!993l) . The uncertainty 
for z, k and a, and the correlation coefficient p z k were used to 
calculate the errors in the magnitudes, Sect. 12.1.41 

In Fig. |2] we show the residuals between our ph otometr y and 
the sta ndard values from ISchuster & Nisser] d!988l) and [Olsen 
(119931) as a function of our calibrated (v-y) colour. We note that 
there are no trends with colour. 



2.1 .3. Photometry of the stars in the science images 

Instead of co-adding the science images, we did aperture pho- 
tometry on each of the images separately. Co-adding the images 
would be difficult since the seeing varied from night to night and 
because of the necessity to apply the extinction correction for 
each night separately. As described later in this section the final 
flux for each star in each filter was calculated using a weighted- 
mean of the individual measurements. This enables us to do a 
more detailed study of the night to night quality. As a quality 
check we compared the flux for the brightest targets for each 
night with a mean flux, calculated for each object for all nights, 
to see if any of the images deviated in flux. None of the images 
deviated. From this we draw the conclusion that the calibration 
was consistent throughout the entire observing run. 

Coordinate lists for the images were created using the task 
DAOFIND in the APPHOT package in IRAF. To establish the 
coordinate lists we used the best y image, since the stars are 
brightest in this filter. We then used this catalogue of coordinates 
for all the other images. We used the aperture photometry task 
PHOT, within the APPHOT package, to measure the flux for all 
objects on the images. 



File name 


Filter 


CCD chip 


Aperture correction 


r556732 


y 


1/2/3/4 


1.265/1.266/1.278/1.266 


r556737 


b 


1/2/3/4 


1.273/1.262/1.274/1.256 


r556748 


V 


1/2/3/4 


1.195/1.191/1.164/1.168 


r556754 


u 


1/2/3/4 


1.171/1.161/1.135/1.148 


r556882 


b 


1/2/3/4 


1.269/1.264/1.259/1.240 


r556887 


y 


1/2/3/4 


1.162/1.166/1.162/1.155 


r556893 


V 


1/2/3/4 


1.268/1.264/1.223/1.238 


r556898 


11 


1/2/3/4 


1.378/1.366/1.312/1.343 


r556899 


11 


1/2/3/4 


1.384/1.362/1.316/1.352 


r556997 


y 


1/2/3/4 


1.254/1.261/1.254/1.244 


r557004 


b 


1/2/3/4 


1.212/1.203/1.204/1.188 


r557009 


V 


1/2/3/4 


1.382/1.382/1.336/1.356 


r557018 


u 


1/2/3/4 


1.209/1.198/1.160/1.185 


r557019 


u 


1/2/3/4 


1.335/1.322/1.282/1.308 



Column 1 lists the file name for the image as named by the observ- 
ing and archiving system on the Isaac Newton Telescope. Column 
2 lists the filter. Column 3 lists the CCD chip number and column 4 
the aperture correction for each CCD in the same order as in column 
3. 



Final magnitudes. Initial magnitudes were calculated for 
each object and night for every image and calibrated for the 
airmass extinction and zeropoint using Eq. (0, but this time 
with subscript s as our calibrated magnitude and, as before, as 
the observed magnitude, with coefficients from Table [3] Since 
the first night did not give us reliable standard star photometry, 
and we thus have no calibration for that night, we normalized 
the magnitudes from that night to the mean of the magnitudes 
for the two following nights. 

For all three nights, erroneous measurements returned from 
PHOT for the individual exposures were removed from the 
data set (i.e. the measurements for which sier, cier and 
pier^O, which are the error in sky fitting, centering algorithm 
and photometry, respectively). The flux was then calculated for 
each star. The final flux, F, was obtained by using a weighted- 
mean flux where the photometric errors returned from PHOT 
(merr) were used as weights. The expression for the final flux is 
thus 



p ^ ; f.jtr) 

* " Z3U Ho) 



(4) 



where the subscript j is the exposure, n is the total number of 
exposures, fj is the flux of the individual exposure and <x, the 
error (merr). These mean fluxes were then converted back to 
magnitudes and the colour terms wer e applied t o get the final 
magnitudes. By definition, y = V (e.g. lOlsenll 19831) . and we will 
henceforth use V instead of y in figures and discussions. 



Aperture correction. When doing photometry on the science 
images we used a fixed aperture of 5 pixels. Applying a curve- 
of-growth we obtained, for each individual image, the aperture 
correction out to 4 x FWHM of the psf. The aperture correc- 
tions were based on measurements of many bright isolated stars, 
typically 20 stars per CCD. The aperture corrections were done 
in flux-space and are listed in Table|4] 



2.1.4. Photometric errors 

The errors in magnitude for the Stromgren photometry were 
calculated using a Monte Carlo Simulation, taking into account 
merr, the uncertainty in zeropoint, extinction coefficient and 
colour term and the strong anti-correlation between the uncer- 
tainty in zeropoint and extinction coefficient. 
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Night 


ky 




Zy 




«i 




Pzk 


14 April 2007 


-0.155 ± 


0.021 


23.009 ± 


0.031 


0.016 ± 


0.005 


-0.92 


15 April 2007 


-0.142 ± 


0.016 


22.989 ± 


0.023 


0.016 ± 


0.005 


-0.96 


Night 


h 




Zb 




a b 




Pzk 


14 April 200/ 


-0.240 ± 


0.031 


23.337 ± 


0.043 


0.009 ± 


0.005 


-0.97 


15 April 2007 


-0.210 ± 


0.036 


23.297 ± 


0.050 


0.009 ± 


0.005 


-0.96 


Night 


kr 




Zr 




a, 




Pzk 


14 April 2007 


-0.349 ± 


0.023 


23.051 ± 


0.034 


0.050 ± 


0.006 


-0.98 


15 April 2007 


-0.353 ± 


0.023 


23.048 ± 


0.032 


0.050 ± 


0.006 


-0.96 


Night 


h 




z„ 




a, 




Pzk 


14 April 2007 


-0.582 ± 


0.036 


23.131 ± 


0.052 


0.065 ± 


0.007 


-0.98 


15 April 2007 


-0.569 ± 


0.013 


23.117 ± 


0.021 


0.065 ± 


0.007 


-0.96 



Columns 2 to 4 list the airmass extinction coefficients, k it zeropoints, z, and colour coefficients, a,, for each filter as indicated with uncertainties 
(Tk, o~z ar, d 0"a respectively. Column 5 lists the correlation coefficient between the uncertainty in extinction coefficient, at, and uncertainty in 
zeropoint, cr.. This coefficient is denoted by p zk . 




Fig. 3. Errors in the photometry for y, b, v and u, as indicated. 
The calculation of the errors is described in Sect. 12.1.41 The few 
stars with anomalously high errors at a given magnitude are also 
discussed there. 

We did this in the following way. For each object, a new 
magnitude was calculated 



m, = niQ i + Am 



(5) 



where mo is the observed magnitude in filter i and Am is a ran- 
dom number, drawn from a normal distribution with a mean of 
and a variance of cr^, where <j m is the error in the magnitude as 
returned from the PHOT task (merr). Additionally, for each ob- 
ject, a new zeropoint and extinction coefficient were calculated 



Zi = zoj + Az 



■Ak 



(6) 
(7) 



where zoj and koj are the zeropoint and extinction coefficient for 
filter i used in the calibration( see Sect. 12.1.2b and Az is a ran- 
dom number, drawn from a normal distribution with a mean of 



and a variance of cr?, see Table [3] Since cr z and o> are strongly 
anti-correlated, we take into account the correlation coefficient 
from the calibration when calculating the random number Ak 
(see Sect. 12.1.21 and Table |3). Finally, a new colour term was 
calculated 



a, = ao i + Aa 



(8) 



where ao,; is the calculated colour term for filter i and Aa is a 
random number, drawn from a normal distribution with a mean 
of and a variance of <x^. 

This process was then repeated 2000 times, generating a new 
set of u, v, b and y magnitudes in each iteration. As the final 
error for each magnitude we adopt the standard deviation of the 
magnitudes from the Monte Carlo simulation. This is calculated 
as follows 



1 2 



(9) 



where the subscript /' is the Monte Carlo iteration number, and 
m the mean magnitude of the distribution. 

Figure [3]presents our final photometric errors, e. The base- 
level error of ~ 0.02 mag (0.035 for the b filter) is due mainly 
to the errors in the photometric calibration, see Table [3] and the 
discussion above. The profile of the errors, i.e. increasing error 
with decreasing magnitude, is dominated by merr (i.e. photon 
statistics). Some stars show a larger e than the majority of stars 
at that magnitude. This is most obvious in the u filter and is a 
statistical feature caused by the number of exposures that are 
included for that star. The u filter suffers most from this since 
stars are generally a lot fainter in this filter, and therefore more 
measurements are rejected due to errors in sky fitting, centering 
and photometry. For example, the stars in u with a larger e for 
a given magnitude have three or fewer individual measurements 
while the stars in the main trend all have five measurements. 
The higher error in b is caused by the higher uncertainty for the 
zeropoint in b (compare Table[3]l. 

2.1.5. Stellar classification using SExtractor 

Contamination by background galaxies is a concern for our 
observations. As our images are uncrowded we can reach 
faint magnitudes with relative ease using aperture photometry. 
However, by using aperture photometry we have no information 
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Fig. 4. SExtractor stellarity index (sclass) as a function of Vo. 
A star takes on values close to 1 and galaxies values close to 
0. The rectangle indicates the area where we investigated each 
object individually (see Sect. 12. 1.5V 



about the shape of the objects. In order to sort the stars from 
background galaxies we used SExtractor v2.5. SExtractor uses 
a tunable neural network trained on realistic simulated images 
to separate galaxies from stars in moderately crowded images 
dBertin & Arnoutsl H996). The SExtractor output, which is of in- 
terest to us, is the stellarity index, sclass; sclass takes on val- 
ues betwe en and 1, where ind icates a galaxy and 1 a star. 
Following [Bertin & Arno uts (1996), objects with a stellarity in- 
dex greater than or equal to 0.5 were identified as stars. 

As can be seen in Fig. |4] SExtactor clearly removes galax- 
ies at bright magnitudes but its ability to distinguish stars from 
galaxies diminishes at fainter magnitudes where the S/N is 
lower. During our membership analysis (see Sect. 14.21 ) we in- 
vestigated objects in a stellarity index range 0.07 < sclass < 0.5 
and magnitude range 1 8 < Vo < 2 1 in order to make sure that 
we did not exclude any objects of interest based on the stellar- 
ity index alone. This area is marked with a solid line in Fig. |4] 
We found two objects of interest, INT 34489 and INT43290 (see 
Sect. |4.2| for a discussion of these two objects). 

2.2. Spectroscopy 

Our spectroscopy was carried out using the multiobject 
spectrograph Fibre Large Array Multi Element Spectrograph 
(FLAMES) at the Very Large Telescope (VLT) on Paranal. 
Operated in Medusa fibre mode, this instrument allows for the 
obser vation of up to 130 targets at the same time (Pasquini et al. 
l2002h . FigureQ]shows the fibre positions on the sky for the stellar 
targets. 21 additional fibres were dedicated to observing blank 
sky. We used the GIRAFFE/L8 grating, which provides a nom- 
inal spectral resolution of R ~ 6500 and a wavelength coverage 
from 821 nm to 940 nm, centred on the Ca n IR triplet lines in 
the spectral region around 860 nm. 



Table 5. Summary of the spectroscopic observations with 
FLAMES. 



Date 


bxp. time [mm] 


15 April 2007 


45 


15 April 2007 


45 


10 May 2007 


45 


17 May 2007 


40 


21 June 2007 


40 


21 June 2007 


40 


21 June 2007 


40 


Total Exp. Time 


295 



Column 1 lists the date of observation and column 2 the exposure 
time. 



2.2.1 . Reduction of spectroscopic observations 

Initially, the FLAMES observations w ere reduced with th e stan- 
dard GIRAFFE pipeline version 2.2 dBlecha et al.ll2000h . This 
version of the pipeline, h owever, was not able to reduce a per- 
sistent glow on the CCD dLind et a l. 2008). This glow then cre- 
ated an extra background that increases towards the red. The ef- 
fect was very large and would have affected the equivalent width 
measurements. Fortunately, during our work with these spectra 
a beta-version of the next version of the pipeline became avail- 
able. The data were thus re-reduced with the GIRAFFE pipeline, 
version 2.5. This pipeline provides bias subtraction, flat field- 
ing, dark-current subtraction, and accurate wavelength calibra- 
tion from a ThAr lamp. It also solved the issue of the CCD glow. 

2.2.2. Spectroscopic measurements 

The 21 sky spectra were combined and subtracted from the ob- 
ject spectra with the task SKYSUB in the SPECRED package in 
IRAF. Figure [5]gives an example of the sky-subtraction process. 

Finally, the object spectra from the individual frames were 
Doppler-shifted to the heliocentric rest frame and median- 
combined into the final one-dimensional spectrum. When com- 
bining the object spectra we used an average sigma clipping al- 
gorithm, rejecting measurements deviating by more than 3 <x, in 
order to remove cosmic rays. 

Radial velocities were determined by a Fourier cross- 
correlation of the combined spectra against a synthetic template 
spectrum using the IRAF task FXCOR. The template consisted 
of three Gaussian absorption lines at the positions of the Ca n IR 
triplet lines, with equivalent widths ( W) representative for red gi- 
ant stars. The radial velocities were determined from a Gaussian 
fit to the strongest correlation peak within a 300 km s~ ! window. 
The uncertainty in the measurement of the radial velocity was 
returned by FXCOR and is based on the Tonry-Davis R-value 
dTonrv & Davislll979l) (see Fig.lTk and b). 

During the Fourier cross-correlation process we performed 
an ocular inspection of the quality of the spectra. Spectra of ob- 
jects fainter than Vo ~ 21.3 showed no clearly visible Ca n IR 
triplet lines and were thus removed from the sample. The S/N 
for these spectra was typically ~ 4 or less. 

The equivalent widths, W, for the Ca n IR triplet lin es were 
measured by fitting a Gaussian profile dCole et al.ll2004l) to each 
of the three lines using the IRAF task SPLOT. From an ocular 
inspection of the spectra we found that the Gaussian profile fitted 
the Ca ii triplet lines better than a Voigt profile. 
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Fig. 5. Example showing the subtraction of sky-emission lines from the stellar spectra. At the top we show a raw, one-dimensional, 
spectrum for one of the brightest stars in our sample, V = 18.4. Below that spectrum is the final spectrum of the sky constructed 
through the combination of all 21 sky fibres. At the bottom we show the stellar spectrum once the sky has been removed. The S/N 
for the final stellar spectrum is 28. The RMS in the sky spectrum is ~ 1.1 counts and ~ 9.8 counts in the stellar spectrum. The three 
Ca n IR triplet lines are indicated by vertical lines in the bottom spectrum. The dashed line indicates the position of zero intensity 
for the final spectrum. 



2.3. Interstellar reddening 

We corrected the photometric magnitudes fo r inters tellar extinc- 
tion using the dust maps by ISchlegel et al.l ( 1998b. Thi s gives 
E(B - V) = 0.062, in agreement with Koch et alJ(l2008bl) . 

IColeman et ail d2007l) used a reddening of E(B - V) = 0.055 
with an uncertainty of 0.005 that represents the variation in red- 
dening over the Large Binocular Telescope 23' x 23' field. In 
Sect. |8.1| we investigate how different values of E(B - V) affect 
the estimated metallicities for the stars. 



We used the Schleg eTet al.l (IT9 98) relations to translate these 
extinction values into the Stromgren system. De-reddened mag- 
nitudes, colours and Stromgren indices will henceforth have the 
subscript 0. 



3. Results 

3. 1 . Colour magnitude diagram in the direction of the 
Hercules dSph galaxy 

Figure[6]presents our colour magnitude diagrams in the direction 
towards the Hercules dSph galaxy. The horizontal branch (HB) 
is seen at Vo ~ 21.2 in Fig. [6^- A large population of foreground 
stars can also be seen with a cut-off at (b - y)o m 0.3, associated 
with the blue limit of the turnoff stars in the Milky Way disk and 
halo. The RGB of the Hercules dSph galaxy cannot easily be 
seen due to the heavy contamination by foreground dwarf stars. 



3.2. Radial velocities 

Figure |7h shows the stellarity index as a function of the derived 
radial velocity for the objects observed with FLAMES. Objects 
with a stellarity index less than 0.5 are excluded from the sam- 
ple (see Sect. 12.1. BY We investigated the objects with a stellarity 
index less than 0.5 and a velocity close to the systemic veloc- 
ity and found that all of these objects are fainter than Vo = 21 
except three objects, INT 43290, INT 40586 and INT 41659. 
INT 43290 lies on top of a background galaxy and is therefore 
contaminated. INT 40586 and INT 41659 have a questionable 
surface distribution on the CCD. We flag them as blends and 
possible double-star systems and identify them in Fig. IT}?. 

Fig-lib shows the magnitude of the stars vs. the derived ra- 
dial velocities. The error- bars on the radial velo cities are from 
the Tonry-Davis R-value (iTonry & Davisl Il979l) estimates. As 
expected, errors in the velocities correlate with magnitude such 
that fainter stars have larger errors. 

In Fig. [8] we compare the distribution of our radial velocities 
in the direction of the He rcules dSph galaxy with the predictions 
of the Besangon model (Ro bin et alj |2003) with the same area 
on the sky as that covered by our observations. The colour and 
magnitude range for the distribution of stars from the Besancon 
model is the same as spanned by the FLAMES targets. As can be 
seen, also in this respect the Hercules dSph galaxy suffers from 
heavy foreground contamination. We note that the observed ve- 
locity distribution is not centred on the field star distribution 
predicted by the model, but is shifted by about -25kms~'. 
Moreover, disregarding the peak created by the Hercules mem- 
bers, the observed distribution is narrower than the model predic- 
tion. This may point to some problem with the model; however, 
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a full investigation of the origin of this difference is beyond the 
scope of the present study. 

4. Finding the giant and horizontal branch stars in 
the Hercules dSph galaxy 

4.1. The ability of the Stromgren c\ index to identify RGB 
stars 

As Figs. [6] and [8] show, the line of sight towards the Hercules 
dSph galaxy is heavily contaminated with foreground stars, mak- 
ing it impossible to determine membership from the colour- 
magnitude diagram alone. Even when radial velocities are added 
the selection remains uncertain, because the mean velocity of the 
Hercules dSph galaxy coincides with the velocity of the (thick) 
disk. 

The c\ index in the Stromgren system gives us the ability to dis- 
entangle the RGB and HB stars in a dSph galaxy from the fore- 
ground dwarf stars. The c\ index is a measure of the Balmer 
discontinuity in a stellar spectrum and is defined as 

ci = (u-v)-(y-b) (10) 

The strength of the Balmer discontinuity depends on the evolu- 
tionary stage of the star. Stars in a plot of ci,o vs. (b - y)o will 
therefore occupy different regions depending on their evolution- 
ary stage. Figure [9] shows which regions are occupied by stars at 
different evolutionar y stages. This classification is adopted from 
ISchuster et alj d2004l) . We also show t he region occup ied by the 
RG B stars in the Draco d Sph galaxy dFaria et alj|2007l) . 

ISchuster et alj ((2004) were mainly concerned with high ve- 
locity dwarf stars and to a lesser extent interested in the redder 
dwarf and RGB stars. As the high-velocity halo stars that they 
studied tend to be fairly blue we will use tracings for dwarf stars 
from Arnadottir et al. (in preparation), and isochrones for RGB 



stars by[Va ndenBer g et alJd20 06) with colour transformations by 
IClem et all d2004h . in the cj.o vs. (b - y)o diagram to define the 
dwarf and giant star regions also in the red. The dwarf sequences 
in Arnadottir et al. (in preparatio n) provide a n extension of the 
preliminary dwarf relation from Olsen (198 4J). The ma jor differ- 
ence between the preliminary relation by lOlsenl d 19841) and the 
new relations is that the new relations are functions of metallic- 
ity. In Fig.|9]we show dwarf sequences for t hree different metal - 
licities. These sequences, in accordance with lOlsenl (fl984). trace 
the lower envelope for the dwarf stars (for (b -y) less than about 
0.55). Note that for redder colours ((b-y) larger than about 0.55) 
all dwarf sequences merge and form a single line that trac es the 
mean values o f the colours. By studying dwarf stars from lOlsenl 
( 19931 fl994llbl) in a ci,o vs. (b - y)o diagram it is possible to 
define an upper envelope for the region occupied by foreground 
dwarf stars. Arnadottir et al. (in preparation) define such an up- 
per envelope. We include this in our plots henceforth. 

To define the RGB region we use two isochrones with 
[Fe/H]= -2.3 and rFe/Hl = 0.0 bvlVandenBerg et all d2006l) and 
colour transformations bv lClem et al.ld2004l) , see Fig.fTOl 

As can be seen from Fig. [10] for giant stars, the c\ index has 
a clear metallicity dependence. This is more prononounced for 
the reddest colours (i.e. the tip of the RGB). However, in spite of 
this, this index still provides a strong discriminant between giant 
and dwarf stars for cooler stars. Note that the dwarf sequence is 
not metallicity dependent at these colours. 

The blue limit for membership determination. Since in the 
c\ x) vs. (b - y)o plane the dwarf and RGB stellar sequences 
converge around (b - y)o ~ 0.5 we need to identify a blue 
limit for stars that we identify as RGB stars. In Fig. [10] we 
show the c\fi vs. (b - y)o for stars in the magnitude range 
15 < Vo < 18. Stars fainter than Vo = 15 are not saturated 
on the images and, based on our membership-determination in 
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Fig. 7. (a) SExtractor stellarity index vs. radial velocity for the objects observed with FLAMES. The error-bars represent the error 
in V ra d as returned by the task FXCOR. The dotted line indicates stellarity index 0.5. (b) Vq vs. radial velocity. • indicates objects 
with a stellarity index greater than 0.5. The error-bars represent the error in V rm { as returned by the task FXCOR. o indicates objects 
with a stellarity index lower than 0.5. Three stars are identified with their INT numbers. These stars are discussed in Sect. l3.2l and 
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Fig. 8. Distributions of radial velocities. The solid histogram 
shows the distribution of velocities from a Besancon model in 
the direction of the Hercules dSph galaxy. The shaded histogram 
shows the distribution of radial velocities for our FLAMES ob- 
servations. Only objects with sclass > 0.5 are included (see 
Sect. 12. 1.5b . Note that the area of the histogram based on the 
Besangon model has been normalized to cover the same area as 
the FLAMES histogram. 



the next section, this magnitude range is bright enough not to 
contain any RGB stars in the Hercules dSph galaxy. We define a 
line that follows the upper envelope of observed dwarf stars to 
separate the RGB stars from the dwarf stars in order to safely 
exclude any dwarf stars. This line is somewhat higher in ci,o 
at a given (b - y)o than the tracing from Arnadottir et al. (in 
preparation). Our selection of RGB stars thus has a blue limit 
that is somewhat colour-dependent. No object bluer than this 
limit will be considered as a RGB star since they have a high 
probability of belonging to the foreground dwarf contamination. 

In Fig.[10]we see 8 stars within the RGB region. A Besancon 
model in the direction of the Hercules dSph galaxy gives 15 
RGB stars in the given magnitude and colour range. The con- 
clusion is therefore that these stars are most likely foreground 
RGB stars belonging to the thick disc of the Milky Way. 

4.2. Membership based on Stromgren photometry. 

We study the ci,o vs. (b — y)o in magnitude bins rather than the 
entire sample of all stars all at once. With this approach the RGB 
region is easier to track. Figure QT| shows the errors in ci_o as a 
function of magnitude. We note the large errors in c\ for fainter 
stars. These errors are due to the higher errors in u at any given 
magnitude. Note that in terms of errors (see Fig. [3), V — 20 cor- 
responds to u ~ 22.5 for our RGB stars. In order to exclude stars 
with large errors, we reject all stars with errors higher than the 
median error, plus 0.015 mag, at a given magnitude (as indicated 
by the solid line in Fig. [TO. We will now use this function in our 
classification of members and potential members of the Hercules 
dSph galaxy. 

Figure [T2h shows a clear and simple separation between the 
foreground dwarf stars and the four stars that fall in the RGB 
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-0.05 and -1 from top to bottom and the dotted line marks the 
upper envelope for dwarf stars (all lines from Arnadottir et al. in 
preparation). 
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Fig. 10. ci o vs. (b - y)o for stars, o, in the direction of the 
Hercules dSph galaxy. Only stars in the magnitude range 15 < 
Vo < 18 are shown. • mark stars t hat fall within the RGB r egion 
as defined by the isochro nes by [V andenBe rg et al] (|2006) with 
colour transformations bv lClem et aiT T2004).~ Metallicities as in- 
dicated. The thin dashed lines indicates dwarf star sequences for 
different metallicities, [Fe/H] = 0.45, -0.05 and -1 top to bot- 
tom and the dotted line marks the upper envelope for dwarf stars 
(all lines from Arnadottir et al. in preparation). The thick line 
marks the empirically determined limit for the foreground con- 
tamination (see Sect. 14. U . 



region. These four RGB stars all have errors in c\ o less than 
the line fitted in Fig.QT] Figure [T2b-d show the three remaining 
magnitude bins. We also note that 9 stars fall within the RHB- 
AGB region, they are discussed in Sect. 14.31 

In Fig. [13] we show e C] vs. Vo and (b - y) for all the stars 
falling in the RGB region. 

In total we find 29 member stars in the RGB region, and 9 in 
the RHB-AGB region. They are listed in Table|9] 



Objects with a stellarity index lower than 0.5 As mentioned 
in Sect. 12.1.51 we found two objects of interest, based on their 
position in the ci o vs. (b - y)o diagram, with a stellarity index 
lower than 0.5. INT 34489 at V = 20.27 has sclass = 0.3 and is 
located right at the edge of the blue limit. An ocular inspection 
of the object indicates that it may be a star, but it is slightly elon- 
gated on the CCD. INT 43290 at V = 19.79 has sclass = 0.35. 
In Sect. |3.2| we found that this object lies on top of a galaxy and 
is therefore contaminated. Given the higher uncertainty of these 
two objects we do not include them in further analysis. The stars 
are listed in Table [6] 

A possible foreground RGB star In Fig.fT2b we found one star, 
INT 42668, that fall within the RGB region but it has a radial 
velocity of V ra d = -2.8 km s -1 (see Sect. [5]l and [Fe/H]c a i = 
-0.5 (see Sect. 18. U . Since it has a radial velocity much lower 
than the mean velocity of the Hercules dSph galaxy, we identify 
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Fig. 11. Errors in the Stromgren colour index cj (see Sect. 14. U . 
The solid line indicates a spline function fitted to fall just above 
the trend with smaller errors. This line is used in Sect. 14.21 and 
Fi g . [T2l and [T3l where we define the RGB stars. 
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Fig. 13. Ec, vs. Vo and (b - y)o for the stars in the RGB region 
identified in Fig. [12] Symbols are the same as in Fig.[T2land the 
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Table 6. Objects with a stellarity index lower than 0.5, a radial 
velocity much lower than the mean velocity of the galaxy or Vo > 
21. 



ID 


RA(2000) 


DEC(2000) 


Vq 




Vrad 


34489 


247.66907 


13.07061 


20.27 ± 


0.02 




43290^ 


247.63172 


12.77969 


19.79 ± 


0.02 


44.9 ±3.1 


42668~ 


247.69964 


12.85175 


19.62 ± 


0.02 


-2.8 ± 1.3 



On top of a galaxy 



Column 1 lists the INT ID. Column 2 and 3 list the coordinates. 
Column 4 lists the Vo magnitude and its associated error. Column 5 
lists the radial velocity and its associated error, both in km s _1 . 



this star as a likely foreground RGB star, thus removing it from 
the sample of member RGB stars. The star is listed in Table[6] 



4.3. Finding the HB and AGB of the Hercules dSph galaxy 

In Fig. [12] we identified RGB members of the Hercules dSph 
galaxy. Additionally, we found 9 RHB-AGB stars. However, it 
is possible to use the Stromgren photometry to further explore 
the RHB-AGB and blue HB stars (BHB). FigureQjJt shows the 
BHB, HB and blue straggler ( BS) regions in t he c\ vs. (b - y) 
plane using the areas defined in Schust er et alJ d2004l) . 

In Fig. [T4b and b we show the previously identified RGB 
and RHB-AGB stars. In addition we show all stars bluewards 
of (b - y)o = 0.5 with an e Cl less than the function in Fig. QT] 
We further divide these stars into two sets according to their V 
magnitude. The HB of the H ercules dSph galaxy i s roughly at 
Vo = 21.2 (compare Fig.loland lColeman et al.l(l2007l) l. We there- 
fore define a box with 21 < Vo < 21.4 and -0.2 < {b -y) < 0.5 
for the potential BHB stars. The second region is given by the 
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Fig. 15. Determination of the distance to the Hercules dSph 
galaxy, see Sect. 14.41 Solid line indicates the median distance 
to the Hercules dSph galaxy. Dashed vertical lines indicate the 
upper and lower sextile. 



magnitude range above and below the first box, i.e. 15 < Vo < 21 
and 21.4 < Vo < 22 with -0.2 < (b - y)o < 0.5. 

Figure [T4h and b then show these stars both in the c\ vs. 
(b - y) and in the colour-magnitude diagram. The stars selected 
in the magnitude range of the HB are all narrowly spaced in 
magnitude and they all fall close to the BHB-RHB sequences. 
Hence we take all 10 of these stars to represent the HB of the 
Hercules dSph galaxy. We find a mean magnitude for the HB of 
Vo = 21.17 ± 0.05 and V = 21.36 ± 0.05. 

4.4. A new distance determination to the Hercules dSph 
galaxy 

We used the relation between the absolute magnitude of the HB 
and metallicity (-2.35 dex a s derived in Se ct. HQ, based on glob- 
ular clusters, as defined in ICarretta et al.l d2000t) to determine 
the absolute magnitude for the HB, Mhb, in the Hercules dSph 
galaxy. Using an apparent magnitude for the Hercules HB of 
'^hb = 21.17 + 0.05 we solve for the distance, d 



«hb - M HB = -5 + 5 • logid) 



(11) 



The errors in the distance determination were calculated us- 
ing a Monte Carlo simulation. For each parameter with an asso- 
ciated error (apparent HB magnitude, metallicity and constants 
in the relation between the absolute magnitude and HB and 
metallicity), a new random apparent HB magnitude, metallicity 
and constants were calculated, and a new distance based on these 
parameters was found. This process was repeated 100000 times. 
In Fig. [15] we show the distribution of the Monte Carlo simula- 
tion. As the final errors for the distance, we adopt the lower and 
the upper sextile of the distribution of simulated distances; the 
lower sextile cuts off the lowest 16.67 per cent of the data and 
the upper sextile cuts off the highest 16.67 per cent of the data. 
We find that Hercules is at a distance of 147+® kpc from us. 
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4.5. Variable stars in the Hercules dSph galaxy 

In Fig. [T4b we see some scatter in the colour-magnitude dia- 
gram above the HB. In agreement with a Besangon model, some 
of these objects are likely HB stars that belong to the Milky 
Way. However, the Besangon model can not explain the over- 
abundance of stars seen at V ~ 20.4 and (b-y)o ~ 0.35. We now 
investigate if these stars are variable stars that, due to variation 
in magnitude, position themselves above the HB. Since we do 
not have more than 3 exposures in the bright filters, y and b, it is 
difficult to determine the variability of these stars as a function 
of time. 

We investigated the variance between the individual expo- 
sures in the y, b, v and u filter, but given the small number of 
exposures we found no strong indication of variability for these 
stars. 

Siegel (2006) report the finding of 15 RR Lyrae variable stars 
in the Bootes dSph galaxy with periods from ~ 0.3 to ~ 0.9 days. 
Since Bootes, with My — -5.8 dBelokurov et al. 2006), is fainter 
than the Hercules dSph galaxy, finding a similar number of RR 
Lyrae stars in our sample is a plausible scenario. 

To further investigate the possible presence of variable stars 
we cross-correlated our photometry with the SDSS photometry 
in order to obtain ugriz photometry for our stars. We then trans- 
formed the ugriz magnitudes on to the Johns on-Cousins UBVR I 
system using the co lour transformations in iJordi et aTl (|2006). 
Bona nos et al.l d2004l) found 163 variable stars in the Draco dSph 
galaxy. Mean magnitudes in V,I and B and amplitudes in V 
for these stars are available on-line. We de-reddened the D raco 
photometry using E(B - V) = 0.027 dBonanos et al.ll2004l) and 
constructe d an ellipse encircling the distribution of the variable 
stars from iBonanos et all (120041) in the Vo vs. (B - V)o colour- 
magnitude diagram (see Fig. [T^). 

Using isochrones by Marig oet al.l (120081) we find that the 
HB for two metal-poor, old populations (-2.3 and -2.0 dex) co- 
incide in V. We therefore normalized the positions of the Draco 
variable stars in Vo so that the ellipse is aligned with the mean 
magnitude of our HB in the Hercules dSph ga laxy, Fig. [T4b . 
We als o include the variable stars in Draco from Bonanos ' et al.l 
(2004) that fall outside the ellipse. We see that the positions of 
the variable stars in the colour-magnitude diagram for the Draco 
dSph galaxy are similar to the positions of our open triangles: 1) 
They have a similar colour range 2) They are more likely to fall 
above the HB than below 3) The spread in Vo at (B - V)o ~ 0.5 
is similar for both data sets. 

Our conclusion from this analysis is that the 8 stars, fainter 
than Vo = 19, likely are variable stars in the Hercules dSph 
galaxy. These stars are listed in Table [7] 

5. Determining the systemic velocity of the 
Hercules dSph galaxy 

Using radial velocity measurements to separate the stars belong- 
ing to a dSph galaxy from foreground stars has proven to be 
an efficient method for membership determination. However, for 
the Hercules dSph galaxy this method is complicated since this 
galaxy has a systemic velocity that falls within the velocity dis- 
tribution of the foreground dwarf stars in the Milky Way. 

In Fig. [8] we show the velocity distribution for our obser- 
vations together with a Besancon model in the direction of the 
Hercules dSph galaxy. As can be seen, the velocity peak of 
the Hercules dSph galaxy lies within the velocity distribution 
of the Milky Way galaxy. A sample of member stars, identified 



Table 7. Possible variable stars in the Hercules dSph galaxy 
identified in Sect. 1431 



ID 


RA(2000) 


rvr > i~> / <-> r\f\f\ \ 

DEC(2000) 


Vo 




(b 


- y)o 


11718 


247.84445 


12.60114 


20.44 ± 


0.02 


0.44 


±0.04 


22960 


247.45022 


12.59089 


20.19 + 


0.02 


0.33 


±0.04 


33388 


247.74828 


12.96788 


20.38 ± 


0.04 


0.13 


±0.06 


41701 


247.78819 


12.78892 


20.74 ± 


0.03 


0.36 


±0.04 


41807 


247.77753 


12.76208 


20.34 ± 


0.02 


0.14 


±0.04 


42113 


247.74934 


12.76745 


20.41 ± 


0.02 


0.25 


±0.04 


42503 


247.71592 


12.77974 


21.78 + 


0.04 


0.05 


±0.06 


43193 


247.64148 


12.79060 


20.54 ± 


0.03 


0.39 


±0.04 



Column 1 lists the INT ID. Column 2 and 3 list their coordinates. 
Column 4 and 5 list the Stromgren magnitude Vo and colour (b-y)o 
and their associated errors, respectively. 



as members based only on radial velocity will therefore contain 
a non-negligible number of foreground stars. However, adding 
knowledge about the evolutionary stage of the stars means that 
we can eliminate the foreground dwarf stars (compare Sect. 14.2b 
and obtain a clean sample which can be used to determine the 
systemic velocity and velocity dispersion for the dSph galaxy. 

In order to illustrate the importance of knowing the evolu- 
tionary stage of the star we will first consider only the radial 
velocities as a means to define a sample of stars belonging to the 
dSph galaxy. After that we will add knowledge about the evolu- 
tionary stage to clean the sample further. 

5. 1 . Using only radial velocities to select RGB members 

We used the ma ximum likelihood method described in 
Wal ker et all |2006) to determine the mean heliocentric veloc- 
ity and internal velocity dispersion for stars in the direction of 
the Hercules dSph galaxy. 

The natur al log arithm of the probability function defined in 
Walker_etaL|(2006) was maximized 

, * , jy , >j ii 
1 v i i 1 v-i v, - u) N 
Hp) = -- X Moj + cr 2 p ) - - V -4 '— - -ln(2n) (12) 

In each iteration we rejected stars deviating by more than 
3 cr as they are likely not members. As can be seen in Fig. [8] 
the objects targeted with FLAMES span a broad range of ve- 
locities. For the first pass through the maximum likelihood it- 
eration we selected stars with lOkms -1 < V ra d < 70kms _1 . 
The maximization converged after 1 iteration. This method gives 
us a mean velocity of 40.87 ± 1.42 km s~' with a dispersion of 
7.33 ± 1.08kms _1 . Stars deviating by less than 3cr from the 
velocity could be considered as possible Hercules dSph galaxy 
members. We find 32 stars in this velocity range. 

5.2. Weeding out foreground dwarf stars with the same 
velocity as the Hercules dSph galaxy 

One of the 32 stars found in Sect. 15 . 1 1 (excluding the three stars 
with sclass < 0.5), INT 42568, lies on the edge of the WFC 
CCD #4 and does not have any photometry available and is 
hence excluded from the following discussions. In Fig. [T6l we 
plot C\fi vs. (b — y)o for the 31 stars considered as possible mem- 
bers based on the radial velocities. We find the following 

- Out of the 31 stars 10 fall on or below the dwarf sequences 
and are therefore excluded. 
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Fig. 16. ci o vs. (b - y)o diagram for stars considered as mem- 
bers from the radial velocity measurement alone (see Sect. 
15. 111. The thin solid lines ind icate the isochrones for RGB stars 
by [V andenBe rget alj d2006l) with colour transformations by 
IClem et all (12004). The thin dashed lines indicate dwarf star se- 
quences for different metallicities, [Fe/H] = 0.45, -0.05 and -1 
top to bottom and the dotted line marks the upper envelope for 
dwarf stars (all lines from Arnadottir et al. in preparation). The 
thick solid line is our lower limit for identification of RGB stars 
(see Sect. [4]). The arrow in the top right corner indicates the mag- 
nitude and direction of the de-reddening applied to the data (see 
Sect. 1231. 
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Fig. 17. Colour-magnitude diagrams, (a) Objects targeted with 
FLAMES are shown as large •. Small dots indicate stars from 
Fig. [6^. (b) Stars that have radial velocities within 3<x of the 
systemic velocity, (c) Stars that are members based on the evo- 
lutionary stage and have measured velocities within 3<x of the 
systemic velocity (see Sect. 15.2b . 



- Of the remaining 21 stars, 18 are RGB members, 2 fall in the 
RHB-AGB region. 

- The last star falls below the RGB solar isochrone but is red- 
der than the blue limit (compare Fig. fT2l and [Toll . This star 
is marked in Fig.[l6](INT 14012). We consider this star as a 
member (but see discussion below). 

Hence, we find about 30 per cent contamination by fore- 
ground stars in our sample. For the 19 RGB stars with the right 
evolutionary stage, we re-derive the mean velocity and disper- 
sion using the maxiumum likelihood method. We note that star 
INT 14012 was rejected during the iteration due to the 3<x limit. 
This star is not included in the final sample of RGB stars. We 
note that one of the remaining 18 stars, INT 42 170, has a ve- 
locity just outside the re-derived more narrow 3<x limit. Since it 
falls short of this limit by only 1.32 km s~' and the error in the 
velocity for the star is 4.77 km s , we keep it in the final sample. 
This star could also be a binary, which would explain its deviat- 
ing velocity. We found no stars that are likely members based on 
the Stromgren photometry, but non-members based on the radial 
velocity measurement. 

We find a mean final systemic velocity of 45.20+ 1 .09 km s _1 
with a dispersion of 3.72 + 0.91 kms -1 . Our conclusion is that 
all stars, except INT 14012, with the right evolutionary stage fall 
within 3<r of the systemic velocity. 

Figure [T7h to c show the colour-magnitude diagrams for ob- 
jects targeted with FLAMES, stars identified as members from 
the radial velocity measurement only and, finally, stars identi- 
fied as members from the radial velocity measurement and using 
photometry to weed out foreground dwarf stars. 

6. Final sample 

Our final sample of RGB, AGB and HB Hercules member stars 
was defined as follows: 1) first we select the RGB stars with the 
right evolutionary stage and a radial velocity within ~ 3<x of the 
systemic velocity (compare Sect. I5.2l and Fig. fTTb). 2) To these 
stars we add stars without spectroscopic measurements, which 
were selected as RGB, AGB, or HB stars as determined from 
photometry (compare Sect. 14.2b . These stars are listed in Table 
[9] Figure [T8h and b show the colour-magnitude diagrams for the 
final sample. 

7. A comparison with previous velocity 
determinations 



Simon & Geha (120071) have obtained radial velocities for 86 
stars in the direction of the Hercules dSph galaxy. In order 
to avoid including foreground dwarf stars they used measure- 
ments of the strength of the Nai lines at 818.3 and 819.5 nm 
to distinguish betwee n dwarf and giant stars. As discussed in 
Schi avon et alj ( 1 1997b the strength of these lines depend on the 
gravity of the star. Hence, it enables a distinction between dwarf 
and giant stars. Out of the 86 stars 29 were identified as members 
based on the measurements of the strength of the Na i lines. 

There are 21 stars in common, i ncluding field dwarf s tars, 
between our spectroscopic study and Simo n~& Gehal (2007H In 
Fig. [19] we show the difference, in units of Gaussia n <x, in stel- 
lar velocity for these stars (see iKlevna et al.ll2002l for a simi- 
lar study of the Draco dSph galaxy). We note that the veloc- 
ity measurements are in agreement, and that there are no ob- 
vious outliers, hence there is no evidence for binaries in this 

2 J. Simon and M. Geha have kindly provided the relevant data so 
that we could do this analysis. 
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Fig. 18. Final colour-magnitude diagrams for the Hercules dSph 
galaxy, (a) Vo vs. (b-y)o- (b) Vo vs. (v— y)o- • are members based 
on the evolutionary stage and have measured velocities within 
3cr of the systemic velocity (see Sect. l5.21 >. o are stars that do not 
have radial velocity measurements but are members according 
to their evolutionary stage (see Sect. 14. 2\ . The filled triangles 
are stars identified as HB stars (see Sect. 14.31) . Her-2 and Her-3 
are discussed in Sect. 18.31 Note that there are two stars in (a) at 
V ~ 18.5 that split in (b). 
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Fig. 19. Difference, in units of Gaussian cr (see Sect. [7), be- 
tween our measured stellar velocities and the velocities from 
ISimon &~G eha (2007). Note that this histogram includes both 
members and non-members. The solid line indicates a Gaussian 
with cr = 1 . 



Fig. 20. Strom g ren ci .o vs. (b - y)o diagram for the stars from 
iimon &Gehal (120071) . • are stars considered as RGB mem- 
bers and o are stars identified as non-RGB members based on 
our Stromgren photometry (see Sect. I4.2l i. x are stars that have 
Vo >2 1 , hence we have not considered their evolutionary stage 
as they are too fa int. The thin solid fines in dicates the isochrones 
for RGB st ars bv|Va ndenBei' g et al.l (120061) with colour transfor- 
mations by IClem et al.l d2004l) . The thin dashed lines indicates 
dwarf star sequences for different metallicities, [Fe/H] = 0.45, 
-0.05 and -1 top to bottom and the dotted line marks the up- 
per envelope for dwarf stars (all lines from Arnadottir et al. in 
preparation). The thick solid line is our lower limit for identifi- 
cation of RGB stars. The arrow in the top right corner indicates 
the magnitude and direction of the de-reddening applied to the 
data (see Sect. 12. 3I >. 



sample. However, given the limited time sampling of the com- 
bined spectroscopy, and the fact that there are stars in our 
spectrosc opic study for which there are no velocity measure- 
ment in ISimon & Gehal d2007l) . we cannot rule out the pres- 
ence of binaries in our study. The presence of binaries could 
inflate the observed v elocity dispersion by a small amount (e.g. 
lQlszewskietailll996l) . 

For the star s in co mmon between our photometric study and 
Simon & Geha (120071) we show a c\ vs. (b - y) diagram in Fig 
l20l From this comparison we find the following 



- Out of the 29 stars considered as members in Sim on & Gehal 
d2007l) . 2 fall in the RHB-AGB region and are therefore ex- 
cluded as RGB members. 

- Of the remaining 27 stars, 12 have Vo > 21, hence we have 
not considered their evolutionary stage as they are too faint. 
4 out of the 12 stars with Vo > 21 fall outside the limits of 

- 2 of the remaining stars fall on or below the dwarf sequences 
and are therefor excluded as members of the Hercules dSph 
galaxy. 

- Thus 13 stars remain that are considered as RGB members 
based on Stromgren photometry. 
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Fig. 21. mi o vs. (v-y)o for the Hercules dSph galaxy. • are mem- 
bers based on the evolutionary stage as derived from photometry 
and they have radial velocities within 3cr of the systemic veloc- 
ity for the Hercules dSph galaxy (see Sect. 15.21 ). The o are stars 
that do not have radial velocity measurements but are members 
according to our photometric criteria as developed in Sect. 14.21 
The solid lines are iso-metallicity lines based on Eq. ( TT4b with 
metallicities as indicated. The shaded area indicates the distri- 
bution of RGB stars from lFaria et alj {2007) for the Draco dSph 
galaxy, adopted from their Fig. 19. 



For the 13 RGB members in common be tween us and 
Simon & Geha d2007l) . using the velocities from Simon & Gehal 
d2007l) . we find a mean systemic velocity of 46 . 1 + 1.30 km s~ 1 
with a dispersion of 4.01 ± 1.08 km s~'. 

For 10 of the 13 RGB members we have determined radial 
velocities (see Sect. 12.2.21 . Given th e small offset in velocit y de- 
termination between our study and ISimon & Gehal (120071) . see 
Fig. [19] it is possible to include the velocities of the 3 stars for 
which we have not obtained a velocity measurement into our cal- 
culation of the mean velocity. We find a mean systemic velocity 
of 44.95 ± 1.02 km s~' with a dispersion of 3.84 ± 0.85 km s~' 
which is in agreement with our result from Sect.P 



There are 15 stars in our final sample of RGB members 
(based on the ev olutionary stage) that d o not have radial veloci- 
ties measured in Simon & Geha] (|2007j). 

Measuring the strength of th e Nai to exclude foreground 
dwarf stars is valid for (V-I) > 1 dGilbert et al.l2006HKoch et al l 
2008c). This corresponds to (b - y) > 0.55 and is confirmed by 
our comparison where the 16 stars that are bluer than this limit 
are identified as either foreground dwarf stars or HB stars by the 
Stromgren indexes. 



8. Metallicities 

8.1. Metallicities based on Stromgren photometry 
8.1 .1 . Determination of metallicities 

The Stromgren filters have proven useful to estimate stella r 
metallicities for RGB stars via the m\ index (e.g. Richtler 1989), 



ni\ — (v — b) — (b — y) 



(13) 



A review of the Stromgren metallicity c alibrations for RGB 
stars, available at the time, is provided by iFaria et alj d2007l) . 
They adopt the lHilkerl (H000) calib ration for thei r RGB stars in 
the Draco dSph galaxy. Since the iHilk er (2000) calibration is 
not valid for [Fe/H] < -2.0 dex and the newly found ultra- 
faint dSph galaxies such as Hercu les are metal-poor, we adopt 
the semi-empirical cali bration by Calamida etail (|2007) onto 
the metallicity scale of IZinn & West) (1 19841 as this calibration 
is valid at least down to -2.4 dex. 

Figure|2T|shows the m^o vs. (v - y)o plane for the stars iden- 
tified as RGB stars in the Hercul es dSph galaxy. 

Equation (fT4l . adopted from Calamida et al.l (120071) . is used 
to convert mi t o and (v - y)o to [Fe/H]c a i for the RGB stars. 



[Fe/H] Cal 



(mi +bi(v-y) + b 2 ) 
(b 3 (v -y) + b 4 ) 



(14) 



where b { = -0.521 + 0.001, b 2 = 0.309, b 3 = 0.159 ± 0.001 and 
b\ — -0.09 + 0.002 (not e that b 2 does not have an error estimate 
in lCalamida et al.ll2007l) . 

For comparison we re-c alculated the stella r metallicities for 
the Draco RG B sample f rom | Faria et all d2007l) using the calibra- 
tion bv lCalamida et alj (120071) (see Fig. |2TT >. We found that the 
Draco dSph galaxy has a me an metallici t y of - 2.0 dex instead 
of -1.74 dex as calculated bv lFaria et al.1 (120071) . A comparison 
between the metallicities derived using the older Iffilkerl d2000l) 
calibration and [Fe/H]c a i shows an offset of ~ 0.3 dex, where 
[Fe/H]cai is more metal -poor. 

In Sect. 12.31 we described how we correct the photometric 
magnitudes for interstellar extinction using E(B - V) = 0.062. 
As a test we re-calculated [Fe/H]c a i using E(B-V) = 0.032 and 
E(B - V) = 0.092, which correspond to an uncertainty of +0.03 
mag. We found that E(B - V) = 0.032 increased and E(B - V) = 
0.092 decreased the metallicity by ~ 0.1 dex and ~ 0.12 dex, 
respectively. 

8.1.2. Error estimates 



Following Faria et al.l d2007l) . the errors in [Fe/H]c a i were calcu- 
lated using a Monte Carlo simulation. 

The simplest version of this approach would be to calculate 
the errors in metallicities taking into account only the measure- 
ment errors from phot PHOT (merr) and propag ate those error s 
through Eq. (fT4b (this is the approach used in lFaria et alj |2007). 
A more involved approach would be to include also the uncer- 
tainties in zeropoints, extinction coefficients and colour terms. 
For comparison we have calculated the errors using both ap- 
proaches. 

In both cases we used the following approach to calculate 
the errors. For every star identified as an RGB member of the 
Hercules dSph galaxy, new random magnitudes (v, b and y) were 
calculated from Gaussian probability distributions with standard 
deviations equal to that of the photometric error for each filter 
using a Box-Muller transformation. m\$ and (v - y)o were re- 
calculated from these new magnitudes and used to derive a new 
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Table 8. Summary of determination of systemic velocities, velocity dispersions and metallicities for the Hercules dSph galaxy. 





Number of stars 


[kms 


-'] 


[kms- 1 ] 


<[Fe/H]> 
dex 


LrC / li. Morli in 
L v l JlVlcQldn 

dex 


This study 
This study 
This study 
S&G 
S&G+d 
Kirby et al. 
Kirby et al.+ci 


28 [RGB stars, ci sel.] 
32 [Only V md sel.] 

18 [RGB stars with V md ] 

29 [Their calculation] 
13 [ci sel. sample] 
20 [Their sample] 

12 [c[ sel. sample] 


40.87 ± 
45.20 ± 
45.0^, ± 
46.1~± 


1.42 
1.09 
1.1 
1.30 


7.33 ± 1.08 
3.72 ± 0.91 
5.1" ±0.9 
4.01 ± 1.08 


-2.35 ±0.31 
-2.34 ±0.30 

-2.58^,+ 0.51 
-2.7~± 0.47 


9 nr+U.I4 
A.iJ_ 03] 

9 9^+0.14 

9 79+O.35 
1 -0.40 



" Value from I Simon & Gehal d2007h 
* Value from lKirbv et al] d2008bl) 



Column 1 and 2 list the number of star s for each study and a short description of how the stars were selected. Note: S&G indicates the sample 
of 29 stars from Simo n & Gehal d2007h . Column 3 and 4 list the systemic velocities and velocity dispersions and their errors. Column 5 lists 
the mean metallicities. Column 6 lists the median metallicities. Errors for the median metallicity are the upper and lower quartile. For the 
calculations we use our metallicites based on Stromgren photometry, but for the Kirbv et al. (2008b) sample we use their metallicities. 




[Fe/H] Cal 

Fig. 23. A comparison between metallicities derived from the m\ 
index and those derived from Ca 11 IR triplet line measurements. 
A[Fe/H] = [Fe/H]cai - [Fe/H] Ca ii. The offset is 0.06 dex with a 
scatter of 0.18 dex. 



metallicity for the star. This process was then repeated 10 5 times 
for each star and since the distributions of the new metallicities 
are not Gaussian around the original [Fe/H], a standard deviation 
based on the sextiles (which is equivalent to lcr in the case of 
a Gaussian distribution) was calculated for each star from the 
distribution of simulated [Fe/H]. This sextile standa rd deviation 
is our error in [Fe/H]c a i (compare iFaria et al] d2007l) ). 

In Fig. [22k we show the errors in [Fe/H]c a i as a function 
of Vo. Note that the five faintest stars have significantly larger 
errors. Also note that the errors are much larger when the uncer- 
tainties in zeropoints, extinction coefficients and colour terms 
are taken into account. We find that it is the errors in the b fil- 
ter that is the largest contributor to the error in [Fe/H]c a i when 
the uncertainties in zeropoints, extinction coefficients and colour 
terms are taken into account (compare TableO. The errors in the 
b filter account for almost 50 per cent of the total error. 

8.2. Metallicities based on Ca 11 IR triplet lines 
8.2.1. Determination of metallicities 

For those RGB stars for which we have FLAMES observa- 
tions we also determined metallicities ([Fe/H]c a n) from mea- 
surements of the equivalent width of the Ca n IR triplet lines 
at A = 849.8, 854.2 and 866.2 nm. We foUow IRutledge et all 



(Il997bl) in defining the line strength of the Ca n IR triplet lines 
as the weighted sum of the W, with lower weights for the weaker 
lines 



z 



W = 0.5 • Wi + W 2 + 0.6 • W 3 



(15) 



where W\ , W2 and W3 are the widths of the individual Ca II IR 
triplet lines in the order of increasing waveleng th. 

As discussed e.g. in IRutledge et alj dl997al) . the strength of 
the Ca n IR triplet lines depend not only on metallicity but also 
on the surface gravity and effective temperature of the star. It is 
possible to remove the effect of gravity and temperature to first 
order by taking into account the position of the star on the RGB. 
This is done by defining the reduced W as 



W 



'=Z 



w + o.64(±o.02)(y-y H B) 



(16) 



where (V - Vhb) is the difference between the V magnitude 
of the star and the V magnitude of the horizontal branch 
(Vhb)- The final rFe/H ] r a ii we re calculated using the calibra- 
tion by Rutledge et al.l dT997al) onto the metallicity scale of 
ICarretta & Grattonl dl997l) . The IRutledge et all dl997al) calibra- 
tion reads as follows 



[Fe/H] CaII = -2.66(±0.08) + 0.42(±0.02)W 



(17) 



Of the 18 RGB stars considered as members based on the 
evolutionary stage and which have measured velocities within 
3cr of the systemic velocity, 3 have too low S/N for measure- 
ments of the W and are therefore excluded from this metallicity 
determination. 

Figure [23] shows a comparison between our two metallic- 
ity estimates. We find that [Fe/H]c a i is on average 0.06 dex 
larger than [Fe/H]c a ii with a scatter of 0.18 dex. This scat- 
ter is consistent with the typical measurement uncertainty of 
[Fe/H]c a ii (see below). In conclusion, the agreement between 
the photometric and spe ctroscopic meta ll icities is very good. 
Since the calibration by ICalamida et al.l (|2007) is valid to at 
least [Fe/H] =-2.4, we conclude that the extrap olation for the 
[Fe/H ] c a 11 to [Fe/H] ~ -2.4 is valid (compare Battaglia et al.l 
120081) . We note that the abundanc e scale of lZinn & Westl d 19841) . 
as used bv lCalamida et al.ld2007l) . in general has ~ 0.2 dex lower 
metallicities at around - 2.0 dex th an the abundance scale of 
ICarretta & Grattonl d 19971) as used bv lRutledge et all d!997al) . 
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Fig. 24. Stromgren ci_o vs. (b — y)o diagram for the stars from 
Kirbv et al.l d200 8b). • are stars considered as RGB members 
and o are stars identified as non-RGB members based on our 
Stromgren photometry (see Sect. 14.21 ). x are stars that have 
Vo >2 1 , hence we have not considered their evolutionary stage 
as they are too faint for our membership determination based on 
Stromgren photom etry. The thin sol i d line s indicates isochrones 
for RGB st ars by | VandenBerg et alj d2006l) with colour transfor- 
mations by IClem et alj d2004l) . The thin dashed lines indicates 
dwarf star sequences for different metallicities, [Fe/H] = 0.45, 
-0.05 and -1 top to bottom and the dotted line marks the upper 
envelope for dwarf stars (all lines from Arnadottir et al. in prepa- 
ration). The thick solid line is our lower limit for the identifica- 
tion of RGB stars. The arrow in the top right corner indicates the 
magnitude and direction of the de-reddening applied to the data 
(see Sect. 12. 31 . 



8.2.2. Error estimates 



Following Battaglia et al.l d2008l) . we define the error due to ran- 
dom noise in the measurement of the W as 



AW 



Vl.5 FWHM 
S/N 



(18) 



where FWHM is the Gaussian full-width-half maximum. The 
S/N varies from ~ 24 for the brightest star to ~ 5 for the faintest 
star. The errors in [Fe/H]c a n were calculated following a Monte 
Carlo simulation procedure similar to the one used in Sect. 18.1.21 
The process was repeated 100000 times. As the final error for 
[Fe/H]caii, we adopt the standard deviation calculated from the 
distribution of simulated [Fe/H]c a ii- In Fig. l22b we show the 
errors in [Fe/H]c a ii as a function of Vq. 



8.3. A comparison with metallicities determined in other 
studies 



Koch et all d2008bl) obtained high resolution spectroscopy of 
two stars in the Hercules dSph galaxy, Her-2 and Her-3. These 
stars correspond to our stars INT 42241 (V = 18.53) and INT 




[Fe/H] Cal 

Fig. 25. A comparison betw een our me t allicitie s derived from 
mi and metallicities from iKirby et all (|2008b). A[Fe/H] = 
[Fe/H]cai - [Fe/H] Kllby . 



41082 ( V = 18.86). The stars are marked in Fig. [18] Koch et al.l 
(12008b!) find [Fe/H] = -2.02 and -2.04 for Her-2 and Her-3, re- 
spectively. 

In Sect. |4.2| we identify both of these stars as RGB stars and 
members of the Hercules dSph galaxy. We only have a radial 
velocity for one of the stars, INT 41082 (Her-3). The velocity of 
this star falls within 3<x of our final systemic velocity (see Sect. 
E3. For this star we find [Fe/H] Ca i = -1-93 and [Fe/H] Ca ii = 
-1.88. For INT 42241 (Her-2) we derive [Fe/H] Ca i = -1.96. 

Hence, there is a difference of 0.11 and 0.06 dex, respec- 
tively, when high resolution spectroscopy Fe abundances and 
metallicities derived from Stromgren photometry are compared. 
This must be regarded as excellent agreement given the com- 
plexities in analysing spectra of such cool, evolved giant stars 
and the general simplifications made when using calibrations of 
photometric measurements to obtain estimates of stellar metal- 
licities. It should also be noted that an overestimate in the red- 
dening of a few hundredths would easily account for this differ- 
enc e (compare Sect. 18 . U . 

IKirby et al. (2008b) studied 20 stars in the direction of the 
Hercules dSph galaxy. Their metallicities are based on a recently 
developed automated spectrum synthesis method tha t takes the 
inform ation in the whole spectrum into account (IKirby et al.l 
2008a). The method was originally developed for globular clus- 
ters in the Mil ky Way and was th en applied t o ultra-faint 
dSph g alaxies in Kirbv et all (l2008bl) . [Fe/H] from IKirby etaf] 
(l2008bl) will henceforth be referred to as [Fe/H]KM>y We have 
cross-correlated our photometry with their 20 stars and found 
the following (see Fig.l24lFI 



- Out of the 20 stars considered as members by IKirby et al.l 
d2008bl) . 2 fall in the RHB-AGB region and are therefore 
excluded as RGB members. 

- Of the remaining 18 stars, 4 have Vo >21, hence we have 
not considered their evolutionary stage but all indications are 
that at least three of them are foreground dwarf stars, com- 
pare Fig. [24] 

- 2 of the 14 remaining stars fall on or below the dwarf se- 
quences and are therefore foreground dwarf stars. 

For the remaining 12 stars we find a mean metallicity of 
[Fe/H] Ca i = -2.25 ± 0.20 de x (Sect. 18711 and a m ean metallic- 
ity based on the values from IKirby et al.l (|2008b) [Fe/H]Kk = 
-2.70 + 0.47 dex. Figure [25] shows the difference between 
[Fe/H]c a i and [Fe/H]icirby We find that the median offset is 



3 E. Kirby has kindly provided the relevant data so that we could do 
this analysis. 
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0.67 dex, where [Fe/H]c a i is more metal-rich. We have inves- 
tigated A[Fe/H] as a function of both magnitude and radial ve- 
locity but found no trend. 

ISimon & Gehal {2007) obtained metallicities from the Ca n 
IR triplet lines and found a mean m etallicity of [Fe/H]=-2.27 
dex. Finally, IColeman et al.l (120071) report a metallicity of 
[Fe/H] =-2.26 dex based on fitting an isochrone to a colour- 
magnitude diagram. 

In summary, we have five different determinations of the 
mean metallicity for the Hercules dSph galaxy. These deter- 
minations range fro m about -2.25 dex (Si mon & Gehal 1 2007; 
IColeman et alj2007l this study) to -2.7 dex dKirbv et al.l20 08b). 
All determinations based on measurements of the Ca n IR 
triplet lines and the mi -index agree on a mean metallicity of 
about -2.3 dex. These determinations also agree with the high- 
res olution analysis by iKoch et al.l (l2008bl) . The determinations 
bv lKirbv etal . (2008b) are ~ 0.5 dex more metal-poor (Table 
|S). This discrepancy is not-negligible. However, f uture high- 
resolu tion studies of additional Hercules stars (as in IKoch et al.l 
(2008b), which essentially is in agreement with our method) 
will determine which methods tend to over- or underestimate 
the [Fe/H]. This has important implications for our understand- 
ing of the origin of the metal poor halo stars (e.g. iHelmi et"aTI 
2006; Koch 200§. 



8.4. Metallicity for the Hercules dSph galaxy 

We find a metallicity range of -2.99 < [Fe/H] CaU < -1.88 with 
a mean metallicity of ([Fe/H]c a i) = -2.35 dex with cr = 0.31 
dex, and a median [Fe/H] of -2.25 dex with an upper and lower 
quartile of -2.1 1 dex and -2.56 dex, respectively. If we exclude 
stars fainter than Vo = 20.5 (compare Fig. l22b) we find a mean 
metallicity of <[Fe/H] Ca i> = -2.26 dex with cr = 0.24 dex. We 
do not detect any significant spatial metallicity gradient in our 
data. 



9. Discussion 

9.1. Spatial distribution 

The new, faint dSph galaxies are in general found to be quite 
elongated (Martin et al. 2008). The dSph galaxy in Hercules is 
no exception, in fact it is one of the new galaxies with the largest 
ellipticity (<? = 0.68). The large ellipticity of these objects might 
be attributable to tidal distortions but could also be due to poor 
sampling statistics. Each galaxy is only repr esented by a limite d 
number of RGB stars in most of the cases (Mai-tineLalJ|2008). 
iMartin et all (120081) and lUral & Wilkinson! (120081) conclude that 
it is entirely possible that the shape determinations for these the 
faintest of galaxies are entirely dominated by shot-noise. This 
does not exclude tidal disruption as an explanation for their, in 
general, very elongated shapes. 

With our Stromgren observations we have searched about 
2/3 of the area on the sky inside the Kin g profile limiting ra- 
dius, as defined by Colema n" et al.l (120071) . for members of the 
Hercules dSph galaxy. Additionally, we have searched about the 
same area on the sky outside the King profile limiting radius 
(see Figs. [Tl and l26li. Through this search we have found in total 
28 stars that are RGB members of the Hercules dSph galaxy. 
These stars are spatially confined to a fairly small, elongated 
area falling inside the core radius from IColeman et all (2007) 
(see Fig.l26h). The RGB members show a slight distortion such 
that there appear to be a few stars trailing the ellipsoid in the di- 




100 



Fig. 27. Distributions of radial velocities. The solid histogram 
shows the distribution of velocities for objects in the direction 
of the Hercules dSph galaxy. The shaded histogram shows the 
distribution of radial velocities for stars identified as RGB mem- 
bers based on the evolutionary stage as derived from photometry. 
The solid- vertical line indicates the systemic velocity as derived 
from the RGB members. 



rection of the Milky Way. These stars are few and none have had 
their radial velocities measured. 

The distribution of the HB, RHB, and potential variable stars 
are shown in Fig.|26*b. These stars confirm the central concentra- 
tion and general shape of the Hercules dSph galaxy. There are 
some HB, RHB and variable stars scattered outside the ellipsoid 
defining the King profile limiting radius. It can not be excluded 
that some of these stars are foreground contaminators. The pres- 
ence of such stars might distort the determination of the shape 
parameters for the galaxy. Measurements of radial velocities for 
these stars would confirm their membership. To our knowledge 
the HB of the Hercules dSph galaxy has not yet been targeted 
for such observations. 

In conclusion we find that for a well-defined sample of 
RGB stars (free from contaminating foreground dwarf stars) the 
Hercules dSph galaxy appears to still have a fairly elongated 
structure, confirming previous studies. 

9.2. Measured velocity dispersion 

As discussed in Sect. [5] the velocity peak of the Hercules dSph 
galaxy lies within the bulk of the velocity distribution of the 
Milky Way galaxy, thus contaminating it. We have shown that 
relying on radial velocites for membership determination yields 
a velocity dispersion of 7.33 + 1.08 kms -1 , whilst excluding 
stars that are not members (i.e. they are foreground dwarf stars) 
yields a velocity dispersion of 3.72 + 0.91 km s _1 . 

It is clear that the velocity dispersion is over-estimated when 
only considering the radial velocity for membership identifica- 
tion. This is expected since the foreground contaminating stars in 
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the direction of the Hercules dSph galaxy have a much broader 
velocity distribution than that for the dSph galaxy (compare Fig. 

In Fig. [27] we show the velocity distribution for our observa- 
tions in the velocity range < V ra d < lOOkms -1 , highlighting 
the objects that are members based on the evolutionary stage as 
derived in Sect. 14.21 As can be seen, objects with a lower veloc- 
ity than the systemic velocity of the Hercules dSph galaxy are 
more likely to be excluded as foreground contaminating dwarf 
stars. Comparing with Fig. [8] verifies that the majority of the ve- 
locities from a Besancon model in the direction of the Hercules 
dSph galaxy have a velocity lower than the systemic velocity for 
dSph galaxy. 

In conclusion we find that when deriving the velocity dis- 
persion for dSph galaxies it is important to have a well defined 
sample of member stars representing the velocity distribution of 
the dSph galaxy. 

9.3. Kinematic sub-structure in Hercules. 



Simon & Geha) d2007l) detected possible evidence of kinematic 
sub-structure in the Hercules dSph galaxy. They found nine stars 
clumped together between 41 and 43 kms~' in a sample of 30 
stars distributed between 30 and 60 kms -1 . 

However, for our final sample of RGB stars in the Hercules 
dSph galaxy a one-sample Kolmogorov-Smirnov test yields a 
significance level of 74% for the null hypothesis that the radial 
velocity distribution is drawn from a Gaussian distribution. The 
conclusion is therefore that we do not see any kinematic sub- 
structure in our sample. 

9.4. Metallicity 



Figure [28a shows the metallicity distribution for the Hercules 
dSph galaxy based on the m\ -index for our final sample of RGB 
members (see Sects. [6] and 18. U . Figure l28b shows the corre- 
sponding normal probability function. The data points only de- 
viate slightly from a linear fit which suggests that the metallic- 
ity distribution is drawn from a single normal distribution. The 
slight deviation from the linear fit is noticeable at low metal- 
licities where the cut-off is less sharp than at high metallicities. 
These features is predicted by models and indicate the occu r- 
rence of strong galactic winds (Lanfranchi & Matteucci 2007). 

In Sect. l8.1.2l we calculated the errors in [Fe/H]c a i- When we 
take into account the uncertainties in zeropoints, extinction co- 
efficients and colour terms we note that the errors are in general 
much larger than the cr of the metallicity distribution (compare 
Fig.l22h where we show the error in [Fe/H]c a i)- However, since 
the uncertainties in zeropoints, extinction coefficients and colour 
terms propagate as a magnitude offset, equal for all stars, we re- 
calculate the errors in metallicities excluding the uncertainties 
in zeropoints, extinction coefficients and colour terms (see Fig. 
l22h). We find that these errors in metallicity, for stars brighter 
than Vq = 20.5, are smaller than the metallicity spread of 0.24 
dex. This enables us to study the star-to-star scatter in [Fe/H]c a i, 
making the profile of the distribution in Fig. [28] significant. We 
conclude that there is an abundance spread in the metallicity dis- 
tribution for the RGB members of up to 1 dex. 



10. Summary 

We have, for the first time, presented a list of Hercules dSph 
galaxy members based on an analysis of radial velocity, evolu- 




[Fe/H] Cal 

Fig. 28. (a) Metallicity histogram for RGB stars in the Hercules 
dSph galaxy. The solid histogram shows the distribution for all 
stars in our final sample. The shaded histogram shows the dis- 
tribution for stars in our final sample with Vo < 20.5 (compare 
Fig. l22k). (b) Corresponding normal probability plot, for stars 
in our final sample with Vo < 20.5, assuming a normal distribu- 
tion. The solid line indicates a linear fit to the data with a mean 
metallicity of [Fe/H]=-2.26 dex and cr = 0.24 dex. 



tionary stage obtained from photometry, and stellar classification 
using SExtractor. In detail we provide the following inventory of 
the Hercules dSph galaxy 

- 28 stars as RGB members based on their evolutionary stage, 
see Sect. I4.2I Of these, 19 have measured radial velocities 
(see Sect. 15.2b . 

- 9 stars as RHB-AGB members based on their evolutionary 
stage (see Sect. 14.2b . Of these, 2 have radial velocities (see 
Sect.ES. 

- 10 stars as BHB members based on their evolutionary stage, 
see Sect.l4~2l 

- 8 stars that are possible variable stars based on their evolu- 
tionary stage, see Sect. 14.51 

Our best determination of the systemic velocity is 45.20 ± 
1.09kms _1 with a dispersion of 3.72 + 0.91 kms -1 . We have 
shown that membership based on radial velocity alone is not a 
good method for the Hercules dSph galaxy, since it has a sys- 
temic velocity that falls well within the velocity distribution of 
the foreground dwarf stars belonging to the Milky Way. 

Stellar metallicities have been determined using the 
Stromgren m\ index with a calibration that translates mi o to 
[Fe/H] for RGB stars. We found a mean metallicity of -2.35 
dex. We also derived metallicities for our stars observed with 
FLAMES from measurements of the equivalent width of the Ca 
ii IR triplet lines. The agreement between the two determinations 
was very good, with an offset of only 0.06 dex. 

Finally, we have estimated the mean magnitude of the HB of 
the Hercules dSph galaxy to Vo = 21.17 + 0.05 based on the 10 
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stars identified as BHB members. This gives a distance of 147^ 
kpc. 
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Fig. 12. c\ o vs. (b - y)o for stars in the direction of the Hercules dSph galaxy. The range of magnitude is indicated at the top of each 
panel. Dots are stars within the given magnitude range. • are stars that fall within the RGB region and have errors less than the line 
defined in Fig.QT] These are members based on their evolutionary stage, o are stars that fall within the RGB and have errors larger 
than the line defined in Fig. QT| The star marked with ® in (d) is flagged as a non-member due to its position in the e Cl vs. (b - y)o 
plane, see Fig. Qj] The star marked with a □ in (c) is a foreground RGB star (see Sect. l4.2l i. The error-bars in (b) represent the error in 
c\ and are only displayed for stars that fall above the line defined in Fig.[TT] The thick error-bars to the right in (c) and (d) represent 
the mean error in c\ for stars that fall above the line defined in Fig. QT|(o). x marks stars that fall in the RHB-AGB region (see 
Sect. 14.31 1. The arrow in the top right corner in each panel indicates the magnitude and direction o f the de-reddening applied to the 
data (see Sect.l23V The solid lines indicate isochrones for RGB stars by VandenBerg et alj d2006l) with colour transformations by 
IClem et al] (120041) . Their metallicities are indicated in (a). The dashed lines indicate dwarf star sequences for different metallicities, 
[Fe/H] = 0.45, -0.05 and -1 from top to bottom and the dotted line marks the upper envelope for dwarf stars (all lines from 
Arnadottir et al. in preparation). The thick line marks the empirically determined limit for the foreground contamination (see Sect. 

S3). 
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Fig. 14. (a) c\$ vs. (b - y)o for the stars in the direction of the Hercules dSph ga laxy. BHB: the blue horizontal branch; HB: the 
remainder of the horizontal branch; BS: the blue stragglers (Schu ster et al.l l2004). The dotted line marks the upper envelope for 
dwarf stars. • are stars identified as RGB members and x are stars identified as RHB-AGB stars based on their evolutionary stage 
(compare Fig.l9land[T2l. A are stars that fall on or are near the HB region but with a large offset to the expected HB magnitude. Filled 
triangles are stars that fall on or near the HB region and with the expected HB magnitude, they are thus identified as HB member 
stars. Error-bars in the top right corner are the mean e c , and e(z,- V ) f° r me stars identified as HB stars, (b) Colour-magnitude diagram 
in the Stromgren system for the stars identified in (a), (c) Colour-magnitude diagram in the Johson-Cousin UB VRI system for the 
stars identified in (a). The solid ellipse outlines the distribution of v ariable stars in the D raco dSph galaxy from Bon anos et al.l 
(2004), normalized in Vq to fit our HB and • marks variable stars from lBonanos et al.l (120041) that fall outside the ellipse. The error 
bar on the ellipse corresponds to the mean amplitude of the variable stars used to define the ellipse. 
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Fig. 22. (a) Errors in the [Fe/H]c a i for all stars identified as RGB members of the Hercules dSph galaxy in Sect. 18.1.21 • indicates the 
errors if only merr and the uncertainties in the coefficients of Eq. ( fT4b are included, o indicates the errors if merr, the uncertainties 
in the coefficients of Eq. ( fT4b . and the uncertainties in zeropoints, extinction coefficients and colour terms are included, (b) Errors 
in the [Fe/H]c a ii as derived in Sect. l8.2.2l for stars identified as RGB members . 
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Fig. 26. Spatial distribution of the stars identified as members from both spectroscopy and photometry. (a)» are members based 
on the evolutionary stage as derived from photometry and they have radial velocities within 3cr of the systemic velocity for the 
Hercules dSph galaxy (see Sect. I5.21 i. The o are stars that do not have radial velocity measurements but are members according to 
our photometric criteria (see Sect. 14.21 ). (b) Filled triangles are HB members. A are probable variable stars and x are stars identified 
as RHB-AGB based on the evolutionary stage as derived from photometry (see Sect. 14.51 and 14.31 respectively). The solid ellipse 
represents the core ra dius and the dotted ellipse the King profile limiting radius of the Hercules dSph gala xy as determ i ned b y 
IColeman et alJ (120071) . Solid lines outline the footprint of the WFC. Central coordinates for the galaxy are from Mar tin et"a"D (|2008). 



Table 9. Hercules dSph galaxy members. 
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Table 9. continued. 



ID RA(2000) DEC(2000) V b v u Chip [Fe/H] Ca i ^[kms' 1 ] 

42621 247.70415 12.80033 20.61 ± 0.03 21.10 ±0.04 21.56 ±0.03 22.50 ±0.05 4 RHB-AGB 

43167 247.64480 12.83356 19.92 ±0.02 20.25 ± 0.04 20.54 ± 0.03 21.35 ±0.03 4 RHB-AGB 

43194 247.64213 12.84056 20.72 ± 0.03 21.09 ±0.04 21.48 ±0.03 22.41 ± 0.05 4 RHB-AGB 



Column 1 lists the ID. Column 2 and 3 list the coordinates. Column 4 to 7 list the Stromgren magnitudes V, b, v and u. Column 8 lists on what CCD chip the star is located. 
Column 9 lists the evolutionary stage of the star. Column 10 lists the photometric metallicity as derived from the Calamida et al.1 (120071) calibration. Column 1 1 list the radial 
velocity, if available. 



